A brief comment on the similarities of the IR solutions for the ghost propagator DSE 

in Landau and Coulomb gauges 

J. Rodriguez-Quintero^ 

' Dpto. Ftsica Aplicada, Fac. Ciencias Experimentales; Universidad de Huelva, 21071 Huelva; Spain. 

This brief note is devoted to reconcile the conclusions from a recent analysis of the IR solutions 
for the ghost propagator Dyson-Schwinger equations in Coulomb gauge with previous studies in 
Landau gauge. 

PACS numbers: 12.38.Aw, 12.38.Lg, ll.15.Tk 



I. INTRODUCTION 



In a series of papers [IHl], we studied the IR behaviour of the solutions for the ghost propagator Dyson-Schwinger 
equation (GPDSE) in Landau gauge. A similar analysis has been very recently carried out in Coulomb gauge [4]. 
This brief note is devoted to reconcile the apparently discrepant pictures resulting from the studies in Landau and 
Coulomb gauges mentioned above. 



II. THE GPDSE IN COULOMB AND LANDAU GAUGES 



As was explained in detail in refs. [IHS]) the low-momentum behavior for the Landau gauge ghost dressing function 
can be inferred from the analysis of the Dyson-Schwinger equation for the ghost propagator (GPDSE) which can be 
written as follows: 
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where F{k'^,^'^) is the ghost propagator dressing function renormalized at the subtraction point, /z^, Hi is the ghost- 
colinear non-perturbative ghost-gluon form factor (usually assumed to be 1) and 
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with a renormalized gluon propagator, 
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where the Schwinger mechanism [7] is invoked to generate, via the fully dressed non-perturbative three-gluon vertex, 
a dynamical gluon mass, M{q'^), that we will approximate by its zero- momentum value in the IR domain ^. 

On the other hand, the authors of ref. [4j recently performed a similar study of the GPDSE in Coulomb gauge. In 
particular, they applied the same strategy followed to investigate the low-momentum Landau-gauge ghost dressing 
function in ref. 1] and took the Gribov's equal-time spatial gluon propagator dressing function 2, 
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^ This is shown to be a good low-momentum approximation for the running mass in ref. [5]. 

^ In very good agreement with the Euchdcan SU{2) lattice results obtained for small lattice couplings in ref. [6]. 
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as the input required to build a kernel and solve the GPDSE, again with the approximation of replacing the fully 
dressed spatial ghost-gluon vertex by the bare one (this is, also in Coulomb gauge, an exact result in the limit of a 
vanishing incoming ghost up to all perturbative orders |5]). Thus, the GPDSE can be rewritten as follows: 
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where / represents the angular integration. 
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It should be emphasized that the ghost propagator dressing function in Coulomb gauge is strictly independent of the 
energy, fcl , as a non-perturbative result of the Slavnov- Taylor identities [9j . 

After assuming a pur e powerlaw behaviour, F{q^) ^ {q^)"'' , for the ghost dressing function and analyzing asymp- 
totically both Eqs. (l]5), one is left in both gauges with the two following well-known cases: (i) ap = {^^ decoupling"^ ) , 
that means zero- momentum finite ghost dressing function; and (ii) ap ^ scaling^'), where the low-momentum 
behavior of the gluon propagator, A{q^) ^ (q^)""^^^ forces the ghost dressing function to diverge at low-momentum 
through the scaling condition: 2ap + (xq — 0. As well in Landau as in Coulomb gauge, a massive gluon propagator 
via the Schwinger mechanism or the Gribov's fomula for the equal-time spatial dressing lead to uq = 1 and thus 
ap = —1/2. In particular, an ellaborated asymptotical analysis of eq. ([I]) leaves us with P] : 
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If ap ^ 0, the perturbative strong coupling defined in the Taylor scheme 
zero- momentum. 
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as can be obtained from Eqs.(3[7|. In the case ap — 0, the subleading correction to the non-zero finite value for the 
zero- momentum ghost dressing function, given by eq. ([t]), is controlled by the well-defined zero-momentum limit of 
Q^tI?^) = {M"^ /q^)aT(q^)-: which is an extension of the non-perturbative effective charge definition from the gluon 
propagator [Tl] to the Taylor ghost-gluon coupling fS^. 

The same two cases result from the analysis of eq. ([5| for the Coulomb gauge in ref. 0] , where a ghost propagator 
dressing function behaving asymptotically as either a constant or F{k^) ^ (A:^)^^/^ is analytically found and confirmed 
by a numerical study. 



III. THE FAMILY OF SOLUTIONS AND THE DIALING PARAMETER 

The GPDSE in eq. ([!]) with the input of a gluon propagator borrowed from lattice QCD calculations is numerically 
solved in ref. [1] and both kinds of solutions in eq. Q were shown to happen controlled by the size of the coupling 
at the renormalization point, g{^i)■ In QCD, one needs to provide a physical scale and a standard manner to proceed 
is by fixing the size of the coupling at a given momentum scale. This can be seen as a boundary condition to solve 
the DSEs. Thus, for any coupling, g{^J,), below some critical value, (?crit, an infinite number of regular or decoupling 
solutions for the ghost dressing, behaving as eq. Q indicates, were found; for g{ii) = ^crit, a unique critical or scaling 
solution behaving as eq. ([t]) was found. It appeared not to be other solutions otherwise. In ref. 1], for a subtraction 
point /i = 1.5 GeV, a critical coupling ^crit — 3.33 and a very good description of ghost propagator lattice data with 
a regular solution of eq. ([I]) for ~ 3.11 were obtained. These results were also recently confirmed [3] by studying 
the coupled system of ghost and gluon propagator DSE in the PT-BFM scheme [T3j. This last work payed attention 
to the critical solution limit by studying how F{0, fj,^) diverges as ^(/x) — >■ (^crit — 1-51, with a subtraction point n = 10 



3 




2 2 

q q 

FIG. 1: (Left) The Taylor coupling defined by eq. ([8| and computed ftom the gluon and ghost propagator results of ref. [3| for 
g^(^i)/(47r) = 0.16,0.17,0.18,0.1812,0.1815,0.1817, with the subtraction point /i = 10 GeV. The curve for the critical limit is 
obtained by applying the results of ref. L, as explained in the text. (Right) The same in a "zoomed" high-momentum domain 
to show all the curves joining each other. 



GeV. One can now apply the perturbative definition of the Taylor strong coupling in eq. (jSj) and compute this coupling 
with the gluon and ghost solutions of [3] in order to see how the critical limit is approached. This is shown in Fig. [T] 
where it can be also seen that all the curves for ax obtained for different values of ^(/i) tend to join each other as 
increases (rigth). As explained in [3], the scaling solution cannot be obtained in the PT-BFM scheme with 
"massive" gluons but only appears as an end-point for the family of regular or decoupling solutions. In Fig. [l] the 
curve for the critical limit is obtained by rescaling, up to giving aT(0) from eq. ([8]) with Hi = 1 at zero-momentum, 
the results at the critical limit for q'^A{q-^)F^{q'^) numerically obtained in ref. [Jf, not by solving the coupled DSE 
system but by applying the lattice gluon propagator as an input to solve the GPDSE. Indeed, the critical value for 
the coupling at /i = 10 GeV can be read from the critical curve in Fig. [I] and one gets (/(/i) ~ 1.56, in fairly good 
agreement with the value of ref. [3]. 

On the other hand, the numerical analysis of eq. (|5| for the Coulomb gauge in ref. |1] also shows both regular 
and critical solution to happen but controlled by i^(0,/i) (or r(0,/i) = 1/F{0,fi)) as a boundary condition with the 
size of the coupling fixed to be g^{fJ.) = g'^ = An x 0.1187 for Nc = 3. In that case, regular or decoupling solutions 
correspond to finite values of i^(G,/i^) and the divergent limit, F{0,fi^) — >■ oo, provides us with the unique critical or 
scaling solution (a similar pattern is claimed to be also found in Landau gauge for the authors of ref. |12|). Again, 
the lattice results for the ghost propagatot seem to agree with a regular solution with F{0, fj?) ~ 10, although much 
larger values for the boundary condition cannot be ruled out. Furthermore, the authors of ref. Q demonstrate that 
the ghost dressing obtained from eq. (jsj) for any F(0,^) behave as c • {k'^)'^^ at asymptotically large k"^ /ni^] where m 
is the Gribov mass scale in eq. (j4]), 7^ is the leading-order ghost-anomalous dimension and c is the same coefficient 
for all the solutions. In other words, in the perturbative domain, the perturbative behaviour is recovered and the 
differences between values of F(fc^,/i^) for different arbitrary inputs of i^(G,/i^) vanish (This can be clearly seen in 
Fig. 2 of [3]). Then, they conclude that the boundary condition is not connected to the renormalization (at least in 
the perturbative regime). 

We are thus left with either a family of Landau-gauge DSE solutions dialed by the size of the coupling at the 
renormalization point or a family of Coulomb-gauge ones dialed by the zero-momentum ghost dressing value as a 
boundary condition not connected to the renormalization. How both pictures can be reconciled? The key point stems 
from the different renormalization prescriptions applied to the ghost propagator in both analyses. In the Landau-gauge 
analysis of refs. [IHS]) the standard MOM prescription, where the Green functions are required to take their tree-level 
expression at the renormalization point and for some particular kinematical choice (this implies F(fj,^,fj?) = 1), is the 
one applied. The prescription applied to the ghost propagator by the authors of ref. [Ij is defined by their eq. (3.20) for 
the renormalization constant Zc{A, [g,r{0)]), where r(0) — 1/F{0, fi^). In particular, this renormalization constant 
depends on the boundary condition, r(0), in such a manner that the value for this boundary condition is rescaling the 
ghost dressing function (and, as can be clearly seen in Fig. 2 of 4J, it does not take the tree-level value, 1, as happens 
in MOM prescription for the subtraction point). In the following, we will show how, depending on the renormalization 
prescription, both patterns can be found for the family of solutions for eq. M in Landau gauge. 
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Let's consider a MOM solution of eq. ([I]) for arbitrary coupling, f;(/i); let's then apply the following transformation: 

g{^i) ^ s g{f,) , Fiq^fi^) ^ -F{q',fi') . (9) 

s 

The properties of eq. ([T]) (the same happens of course for eq. ([5|) guarantee that, for any s being a c-number, the 
transformed dressing function verifies the DSE equation with the transformed coupling (of course, MOM prescription 
implies s = 1). Then, if one chooses s = g/gin) and apply the transformation to every solution of the MOM family, 
we will be left with a one-to-one correspondence between these solutions and the new ones 

F{q\^,^) ^ ^Mpi^q^^,^), (10) 
9 

for the fixed coupling 5, which can be identified by the zero-momentum value, -F(0, fi^). This new family of transformed 
solutions obeys the same pattern of the Coulomb gauge family in ref . |3] . It is interesting to notice that the strong 
coupling defined in the Taylor scheme can be also obtained from the transformed solutions as 

aT{q^) ^ f q^^{q^^,■')F{q\^,^) ^ ^ q^ ^{q\ ^,^)F{q^ ^J?) , (11) 
47r 47r 

alth oug h it is obvious that neither F nor the coupling are in MOM scheme. More interesting enough is to realize from 



eq. ( 11 1 that: (i) F does not depend on /i, as far as one applies the same renormalized gluon propagator to obtain any 
solution for arbitrary coupling, (?(/i), as done in [J]; (ii) the transformed ghost dressing function for the critical MOM 
solution [gip) — > 5crit) corresponds to the scaling solution for the critical boundary condition, F{0,fi'^) 00 and 
diverges as F{q-^, /i^) ~ (9^)"^^^ for q^ — )■ 0, as it clearly results from eq. ([s]); (iii) finally, as the difference between the 
Taylor strong couplings computed for any two arbitrary values of (/(/i) vanishes at asymptotically large momentum 
(see Fig. [1]), the same should happen to F for different values of F(0,/i^). All this is claimed by the authors of [T to 
identify the family of solutions for the Coulomb-gauge eq. 



IV. CONCLUSIONS 



We thus conclude that the behaviour of the family of Coulomb-gauge GPDSE solutions in ref. [4] is analogous to 
the one described in refs. [1-31 for Landau gauge, although not renormalized in MOM scheme but after applying the 
transformation of eq. Q. The input parameter for the solutions in ref. T is the zero- momentum ghost dressing, 
which can be interpreted as a boundary condition and put in connection with the Gribov problem. On the other 
hand, for Landau gauge and MOM scheme, g{^) is related to the strong coupling in Taylor scheme while the size of 
the fixed coupling, g, after applying eq. ^ is physically meaningless. Thus, as done in ref. [3], the critical value for 
g{lJ,) can be used to derive a critical value for which can be compared with lattice evaluations or experimental 
determinations to investigate whether the critical solution can be ruled out. 
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